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Abstract. We present new parameterizations of vector and axial nucleon form factors. We maintain an 
excellent descriptions of the form factors at low momentum transfers, where the spatial structure of the 
nucleon is important, and use the Nachtman scaling variable £ to relate elastic and inelastic form factors 
and impose quark-hadron duality constraints at high momentum transfers where the quark structure 
dominates. We use the new vector form factors to re-extract updated values of the axial form factor 
from neutrino experiments on deuterium. We obtain an updated world average value from i/ M d and pion 
electroproduction experiments of Ma = 1.014±0.014 GeV/c 2 . Our parameterizations are useful in modeling 
neutrino interactions at low energies (e.g. for neutrino oscillations experiments). The predictions for high 
momentum transfers can be tested in the next generation electron and neutrino scattering experiments. 

PACS. 13.40.Gp Electromagnetic form factors - 13.15.+g Neutrino interactions - 13.85.Dz Elastic scat- 
tering - 14.20.Dh Protons and neutrons - 25.30.Bf Elastic electron scattering - 25.30.Pt Neutrino 
scattering 



1 Introduction 

The nucleon vector and axial elastic form factors have 
been measured for more than 50 years in e~ N and vN 
scattering. At low Q 2 , a reasonable description of the pro- 
ton and neutron elastic form factors is given by the dipole 
approximation. The dipole approximation is a lowest-order 
attempt to incorporate the non-zero size of the proton into 
the form factors. The approximation assumes that the pro- 
ton has a simple exponential spatial charge distribution, 
p(r) — poe~ r / r °, where r$ is the scale of the proton radius. 
Since the form factors are related in the non-relativistic 
limit to the Fourier transform of the charge and magnetic 
moment distribution, the above p(r) yields the dipole form 

defined by: G^ A (Q 2 ) = C v > A /(l + jjjr^f ■ Here C V > A = 

(l,g A ), 9A = -1.267, My = 0.71 (GeV/c) 2 , and M A is the 
axial mass. 

Since Ma is not equal to My, the distribution of elec- 
tric and axial charge are different. However, the mag- 
netic moment distributions were assumed to have the same 
spatial dependence as the charge distribution (i.e., form 
factor scaling). Recent measurements from Jefferson Lab 
show that the ratio of ^ Ep falls at high Q 2 challeng- 

<JMp o -v o 

ing the validity of form factor scaling [T] and resulting in 
new updated parameterizations of the form factors [2], [3]. 
In this paper we present parameterizations that simulta- 
neously satisfy constraints at low Q 2 where the spatial 
structure of the nucleon is important, and at high Q 2 



where the quark structure is important. A violation of 
form-factor scaling is expected from quark-hadron duality. 
We use our new vector form factors to re-extract updated 
values of the axial form factor from a re-analysis of previ- 
ous neutrino scattering data on deuterium and present a 
new parameterization for the axial form factor within the 
framework of quark-hadron duality. 



2 New Parametrization 

The new parameterizations presented in this paper are 
referred to as the duality based "BBBA07" parameter- 
ization. Our updated parameterizations feature the fol- 
lowing: (1) Improved functional form that adds an addi- 
tional Q 2 dependence using the Nachtman scaling vari- 
able £ to relate elastic and inelastic form factors. For elas- 
tic scattering (x = 1) £P> n ' N = 2 where 

r p , n ,N = Q 2 /4Mp nJv . Here M p ^ N are the proton (0.9383 
GeV/c 2 ), neutron (0.9396 GeV/c 2 ), and average nucleon 
mass (for proton, neutron, and axial form factors, respec- 
tively). (2) Yield the same values as Arrington and Sick 
[4J for Q 2 < 0.64:(GeV/c) 2 , while satisfying quark-hadron 
duality constraints at high-Q 2 . 

For vector form factors our fit functions are Ajv(£) 
(i.e. A Ep (e), A Mp (?), A En (C), Am nip 1 )) multiplying 
an updated Kelly [3] type parameterization of one of the 
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Fig. 1. Ratios of Ge p (a), Gm p /^ p (b), Geu (c) and GMn/i-in (d) to Gd- The short-dashed line in each plot is the old Kelly 
parameterizations (old Galster for Gb„). The solid line is our new BBBA0725 parameterization for - = 0.0, and the long-dashed 
line is BBBAO743 for - = 0.2. The values of £ and the corresponding values of Q 2 are shown on the bottom and top axis. 
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Table 1. Parameters for G B ^ v and G M ^ y upd 



Our parame- 
terization employs the as-published Kelly parameterization to 
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Ep and an updated set of parameters for G MP 
that includes the recent BLAST 6 results. 



proton form factors. The Kelly parameterization is: 

Em u 
k=0 a kT p 



(Q ) 



G Kelly (Q 2 ) 



where do = 1 and m = 1 

In our analysis, we use all the datasets used by Kelly [3], 
updated to include the recent BLAST [6] results, to fit 
Ge p , Geu, Gmp/^v an d Gmu/^u (Mp = 2.7928, \i n = 
— 1.9130). Our parameterization employs the published 

Kelly functional form to G^ ly , and an updated set of 
parameters for G*, e p y ~ up (Q 2 ). The parameters used for 



G*Ep ly and G^ { e p Uy upd are listed in Table [lj and A N (£) is 



T Ep 

given by 



n 

pj(o=Pi n 



Each Pj is a LaGrange polynomial in £. The £j are equidis- 
tant "nodes" on an interval [0,1], and pj are the fit param- 
eters that have an additional property Apf(£j) — Pj- The 
functional form An(£) (for Ge p , Gm p , Geu, and Gmu) is 
used with seven pj parameters at £j=Q, 1/6, 1/3, 1/2, 2/3, 
5/6, and 1.0. In the fitting procedure described below, the 
parameters of An(£ ) are constrained to give the same vec- 
tor form factors as the recent low Q 2 fit of Arrington and 
Sick [4] for Q 2 < Q.M(GeV/c) 2 (as that analysis includes 
coulombs corrections which modify Ge p , and two pho- 
ton exchange corrections which modify Gm p and Gmu)- 
Since the published form factor data do not have these 
corrections, this constraint is implemented by including 
additional "fake" data points for Q 2 < 0M(GeV/c) 2 . 
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Table 2. Fit parameters for yljv(£), the LaGrange portion of the new parameterization. Note A 2 £ In , A 2 En , and A F 5 A are con- 
strained to have - = at £ = 1, and Aj^ n , j4^ n , are constrained to have - = 0.2. 



Our fits to the form factors are: 

-)2\ /,, a _ _ I CV\ w r~<Ketty—upd/ / ~ l 2\ 



G Ep {Q 2 ) = A Ep (e) x G K E f\Q 2 ) 
Gm„(Q 2 )/m« = A^fiC) x G Mp (Q 2 )/v P 



G En (Q 2 ) = A 2 E 5 f(C) x G Ep (Q 2 ) 



l + br n 



where we use our updated parameters in the Kelly param- 
eterizations. For G En the parameters a=1.7 and b=3.3 are 
the same as in the GalsterJS] parametrization and ensure 
that dG En /dQ 2 at for Q 2 = is in agreement with mea- 
surements. For convenience, we also provide fits for the 
form factors G Ep and Gm p I ' Hp that give very close to the 
same values, but use the dipole form instead: 

G Ep (Q 2 ) = A Ep ^di P oie{S, p ) x GY)(Q 2 ) 

Gm P {Q 2 )/ Hp — AMp-dipole(£ P ) x GY)(Q 2 ) 



G En , 
(Q 2 - 



The values A(£)=pi at £i=0 (Q 2 = 0) for G Mp , G Ep , 



Gmu are set to to 1.0. The value A(£) = 
oo ) for Gm p and G Ep is set to 1.0. 



Z P7 at £,-=1 



The value A(£)=pj at £j=l for Gmu and G En are fixed 
by constraints from quark-hadron duality. Quark-hadron 
duality implies that the ratio of neutron and proton mag- 
netic form factors should be the same as the ratio of the 
corresponding inelastic structure functions in the £=1 

F2p 



limit. (Here Fi = £J^, 



(0) 



r 2 
r 2 



2l> 



1 + 4^ 

U 

4+ £ 



A 2 Mn (t = 1) 



We ran fits with two different values of - at the £=1 
limit: ^ =0 and 0.2 (corresponding to j^ 2 - — 0.25 and 
0.4286). The fit utilizing ± = is Af In , and the fit uti- 
lizing ^ = 0.2 is A^ n . The final parameters for both 
cases of ^ , are given in Table [2] (or download computer 
code|21j). The difference between these two sets is in- 
dicative of the theoretical error of our parameterization. 
Our parameterizations are within the error band of re- 
cent theoretical fits based in dispersion relations;!^. Since 
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Fig. 2. The constraint used in fitting Geu stipulates that 
G 2 En /G 2 Mn — G E p/G 2 Mp at high £. The solid line is ^t E ^ and 

}~ Ep , , and the short-dashed line is 5? E " and !^ E " \ . 



our fits are constrained to give the same vector form fac- 
tors as the recent low Q 2 fit of Arrington and Sick |4] for 
Q 2 < 0.64(GeV/c) 2 , they are in agreement with the ex- 
perimental measurements of the proton and neutron rms 
radii. (Note that as discussed in reference [ID], the nucleon 
rms radius should be determined from fitting a polyno- 
mial of second order to the low Q 2 form factors. The com- 
monly used polynomial of first order yields radius values 
which are too small). 

The value A(£)=pj at £j=l for G En is set by another 
duality-motivated constraint. R is defined as the ratio of 
deep-inelastic longitudinal and transverse structure func- 
tions. For inelastic scattering, as Q 2 — > 00, R n — R p . If 
we assume quark-hadron duality, the same should be true 
for the elastic form factors at £=1 ( Q 2 — > 00) limit: 



R n (x = l;Q 2 
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^ E n 
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Table 3. Ma (GeV/c 2 ) values published by ^ M -deuterium experiments [13] and updated corrections AMa when re-extracted 
with updated BBBA200725 form factors, and g a =-1.267. Also shown is updated Ma from i/jjHydrogcn — >n~ n [14] . 



G2 I (~*2 r^l //^2 

En/ u Mn ~ { - z Ep/ u Mp 

In order to constrain the fits to Geu at high Q 2 we 

G 2 

have assumed that the values of arc the same as the 

G 2 

measured for the three highest Q 2 data points for 

Ge p , and included these three "fake" data points in the 
Geu fits. In addition, the R n — R p condition yields the 
following constraint at £ = 1: 

where b/a — 1.7/3.3. As there are two parameter sets 
^Mn 3 (0; we have produced two parameter sets A 2 ^ 3 as 
shown in Table [2] 

The new form factors Ge p , Gm p /^ p , Gmu/^u, and 
Geu are plotted in Figure 1 as ratios to the dipole form 
Gl. 

As seen in Table v2\ -Ajv(£) is not needed for Gm p as it 
is very close to 1.0. For Ge p it yields a correction of 1% 
at low Q 2 (because it is required to agree with the fits of 
Arrington and Sick 4 (which include two photon exchange 
and Coulomb corrections). For Geu and Gmu it is used 
to impose quark-hadron duality asymptotic constraints. 
Figure |2j shows plots of the data and fits to yg^j and 

rS^r (for the i = at £ = 1 case). 

3 Re-extraction of Axial Form Factor 

Using our updated BBBA20 0725 form factors and an up- 
dated value gA — -1.267, we perform a complete reanalysis 
of published v quasielastic [13] (QE) data on deuterium 
{v^ n — * /i~p) using the procedure described in detail in 
ref. [TT1I2"]. We extract new values of Ma with updated 
form factors (FF) and also include radiative corrections [5] 
(RC). Although of lower statistical significance, for com- 
pleteness we also include all available antineutrino data 
on hydrogen targets [T4] . 



The average of the corrected measurements of Ma 
from Tablets Mf utermm = 1.016 ±0.026 GeV/c 2 . This 
is in agreement the average value of M A ° n =l. 014 ± 0.016 
GeV I c 2 extracted from pion electroproduction experiments 
after corrections for hadronic effects. [T2]. The average of 
the v„ and electroproduction values is 

M world-uverage =l {)U ± Q QU GeV/( ,2^ 

This precise Ma is smaller than the recent results (for 
Q 2 > 0.25 (GeV/c) 2 ) reported by MiniBoonepES] on a car- 
bon target (M A arbon = 1.25 ± 0.12 GeV/c 2 ) and by the 
K2K^9\ collaboration on oxygen (M A V9en = 1.20 ± 0.12 
GeV/c 2 ). Both experiments use updated vector form fac- 
tors. Although the collaborations attribute the larger Ma 
to nuclear effects, there are theoretical arguments that 
Ma in nuclear targets should be smaller [20 than (or the 
same[16 ) as in deuterium. This Ma discrepancy is im- 
portant for v oscillations experiments since it affects the 
normalization (at high energies the QE cross section is ap- 
proximately proportional to Ma) and non-linearity of the 
QE cross section, which is relevant to the extraction of v 
mass difference and mixing angle. 

For deep- inelastic scattering, the vector and axial parts 
of i<2 are equal. Local quark-hadron duality at large Q 2 
implies that the axial and vector components of F% tlc 
arc also equal, which yields: 

[FA(Q 2 )A2=V2? = 

{G V E ) 2 (Q 2 ) ± t n (G v m (Q 2 )) 2 /(\ ± t n ), where 

G V E (Q 2 ) = G Ep (Q 2 ) - G En (Q 2 ), and 

G\[{Q 2 ) = Gm p (Q 2 ) — Gmu(Q 2 )- 

We extract values of Fa(Q 2 ) from the differential cross 
sections using the procedure of ref. [TT] . The overall nor- 
malization is set by the theoretical QE cross section [22 . 
We then do a duality based fit to Fa(Q 2 ) (including pion 
electroproduction data) of the form: 

F a \q 2 ) = Af A (£ N ) x G£(Q 2 ). We impose the con- 
straint A"p A (£i = 0) = p\ = 1.0. We also constrain the fit 
by requiring that Af A (£ N ) yield F A (Q 2 ) = F A (Q 2 ) A2=V2 
by including additional "fake" data points) for £ > 0.9 
(Q 2 > 7.2(GeV/c) 2 ). 
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Fig. 3. (a) Fa(Q 2 ) re-extracted from neutrino-deuterium data divided by Go(Q 2 )[22]. (b) Fa(Q 2 ) from pion electroproduction 



divided by G£(Q 2 



corrected for for hadronic effects [12]. Solid line - duality based fit; Short-dashed line - Fa(Q 2 )a2=\ 



Dashed-dot line - constituent quark model [17] 



Figure 3(a) shows Fa(Q 2 ) extracted from neutrino- 
deuterium experiments divided by G^,(Q 2 )[22]. Figure 3(b) 
shows Fa{Q 2 ) extracted from pion electroproduction ex- 
periments divided by G^(Q 2 )[22]. These pion electropro- 
duction values can be directly compared to the neutrino 
results because they are multiplied by a factor Fa(Q 2 , Ma = 
1.014 GeV/c 2 )/F A (Q 2 , M A = 1.069 GeV/c 2 ) to correct for 
AM a — 0.055 GeV/c 2 originating from hadronic effects [12 . 
The solid line is our duality based fit. The short-dashed 
line is Fa(Q 2 )a2=V2 - The dashed-dot line is a constituent- 
quark model[T7] prediction. 



4 Conclusion 

In conclusion, our new parameterizations of vector and 
axial nucleon form factors use quark-hadron duality con- 
straints at high momentum transfers, and maintain a very 
good descriptions of the form factors at low momentum 
transfers. Our new parameterizations are useful in model- 
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